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Abstract 

The  dielectric  properties  of  soils  and  plastic  targets  are  given.  The  approach  for  finding  exact  pole 
locations  fi'om  a  dielectric  infinite  slab,  a  dielectric  sphere,  and  a  dielectric  infinite  cylinder  is 
presented  (and  applied).  The  exact  poles  for  the  slab  and  sphere  are  found  for  various  sizes  and 
dielectric  properties.  There  is  a  discussion  (and  use)  of  Carl  E.  Baum’s  perturbation  formulas  for 
dielectric  slabs  and  spheres.  George  Hanson’s  results  (formulas,  from  using  Carl’s  same  approach) 
are  given  (and  used)  for  an  infinite  dielectric  cylinder.  The  poles  found  fi'om  the  perturbation 
formulas  are  compared  with  the  exact  poles  and  an  “effective”  region  for  using  the  perturbation 
formulas  is  established. 
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L  Introduction 

The  problem  of  finding  buried  UXO  (unexploded  ordnance)  and  dielectric  mines  has  been  well 
documented  in  recent  years.  According  to  the  Department  of  Defense  there  are  approximately 
nine  million  acres  of  past  and  present  test  ranges  that  contain  ordnance  [MMWG,  1996],  As  a  US 
citizen  what  is  perhaps  most  disturbing  is  that  there  are  between  1200  and  1700  Fonnerly  Used 
Defense  Sites  (FUDS,  now  public  or  private  land)  which  are  believed  to  be  contaminated 
[MMWG,  1996],  As  for  dielectric  mines,  war  zone  areas  like  Bosnia,  Iraq,  and  Afghanistan  have 
demonstrated  just  how  lacking  in  technology  we  are  in  finding  buried  plastic  targets. 

As  far  as  R&D  funding  goes  we  have  reached  a  point  where  the  government  w^ts  something 
that  can  tell  the  user  whether  or  not  the  underground  object  that  is  being  detected  is  in  fact 
ordnance  (and  perhaps  what  type  of  ordnance  it  is)  as  opposed  to  just  being  an  old  soup  can, 
plastic  toy,  or  ordnance  that  has  exploded  and  is  now  environmentaUy  inert.  Well  if  these  things 
were  in  the  air  we  might  be  able  to  do  this;  however,  distinguishing  objects  underground  can  be  an 
electromagnetic  nightmare. 

A  common  cry  fi-om  the  scientist  and  military  personnel  working  on  this  project  has  been  to 
establish  a  library  of  what  the  electromagnetic  pole  signatures  of  the  various  ordinance  and  mines 
look  like.  With  such  a  library  a  user  could  know  if  what  has  been  detected  is  of  interest. 
Developing  this  library,  however,  is  no  simple  task.  In  1994  Carl  E.  Baum  developed  an  approach 
for  using  perturbation  formulas  (much  more  condensed,  yet  approximate,  mathematical  formulas) 
to  find  the  signatures  of  basic  dielectric  geometrical  shapes  in  a  dielectric  medium  [Baun^  1994]. 

It  is  the  purpose  of  this  report  to  establish  the  effectiveness  of  these  formulas  by  comparing  them 
to  the  exact  cases.  If  their  effectiveness  is  proven  then  the  same  approach  could  be  used  in 
creating  our  library. 

n.  Dielectric  Properties  of  Soils  and  the  Ordnance 

Scattering  of  electromagnetic  waves  occur  when  there  is  a  difference  between  the  dielectric 
properties  of  the  target  and  the  medium  [McCann  et  al.,  1988].  Dielectric  targets  of  interest  in  the 
ordnance  and  mine  problem  typically  have  a  permittivity,  e,  of  about  2.5  e^.  The  earth,  _ 

unfortunately,  is  not  as  predictable. 

The  dielectric  properties  of  soil  not  only  depend  upon  the  type  of  soil  but  also  the  amount  of 
moisture  (and  minerals)  in  the  soil.  For  clays  as  the  water  content  rises  fi'om  0  to  40  /o  s  rises 
quickly  fi'om  3  to  [Wang,  1978].  For  sandy  soils  e  rises  fi'om  3  to  [Wang,  1978]. 

The  conductivity,  o,  of  the  soUs  rises  fi'om  =.001  to  =.  1  S/m  for  a  rise  in  water  content. 

For  the  results  found  in  the  follovring  pages  I  will  use  a  target  permittivity,  ,  of  2.5  and 
will  look  at  soil  permittivities,  ,  ranging  fi'om  4  to  20 and  conductivity’s,  a ^ ,  ranging  fi'om 
.001  to  1  S/m.  I  \>rill  also  look  at  the  case  where  the  target  is  in  water  =  81  (and  salt  water 

=  81  cr,  =  4  S/m),  as  well  as  the  case  when  f  i»f2  • 
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nL  The  Exact  Approach 

The  signature  differences  of  the  scattered  waves  from  a  target  not  only  depend  upon  these 
dielectric  differences  but  also  upon  the  various  shapes  and  sizes  of  the  targets.  The  scattered 
wave-forms  for  basic  geometric  shapes  have  been  derived  from  many  different  physicist.  The 
poles  of  these  waves  can  be  found  by  solving  for  the  poles  of  the  scattering  coefficients.  In  the 
following  section  I  will  present  these  derived  wave-forms  and  there  coefficients  for  a  dielectric 
sphere,  infinite  slab,  and  infinite  cylinder. 

Before  I  present  these  three  cases,  however,  I  should  define  a  few  important  terms.  The 
relative  dielectric  constant,  ,  will  be  defined  as 

_  target  _ 
medium  e. 

The  propagation  constants  for  the  medium,  y , ,  and  the  target,  y  2,^^ 


se^ 


1 

2 


and 

r2 


where  5  =  Q  +  j(i)=  the  complex  frequency. 
Finally,  the  ratios  of  these  propagation  constants,  will  be  defined  as 


Ll 
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=  e.  2 
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Sphere: 

The  series  solution  for  a  plane  wave  falling  upon  a  sphere  was  first  introduced  by  Mie  in  1908. 
In  1941  Julius  Adams  Stratton  provided  a  more  extensive  and  complete  model  of  the  problem  in 
his  infamous  book  entitled  Electromagnetic  Theory.  For  derivations  of  the  following  equations  I 
refer  you  to  Stratton,  1941  or  the  Laplace  form  in  Baum,  1992. 

The  scattered  electric,  E,  and  magnetic,  H,  waves  from  a  dielectric  sphere  are 


“  2«  +  l 


E^=EJ(s)1^ 


Y 1  n=\ 
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Where  r  is  the  radius  of  the  sphere  and  a„  and  are  the  scattering  coefficients.  These  scattering 
coefficients  are  functions  of;'  ,r  and  y^r,  and  they  depend  upon  Bessel  functions.  Where  here  I 
will  be  using  the  modified  spherical  Bessel  functions: 

=  +  (-*)”' *.w] 


The  scattering  coefficients  have  been  derived  as 

-z[r2rJn(r2r)]  J„(r^r)-r2rl„(y,r)^ - — 


f 

lrxri„(r,r)] 


-[r2"«0'2^)]  K(rxr)-r2rtAr2r)^^—-y^ — 


r  ..  .l'.,  X  x[rA(r/)] 

Ir2ri„ir2r)]  hiYxr)-Y2rKiy2rr  - 


b.  = 


r  .,  o',,  S  ./  ^[ri^^„(r/)] 

[y2^niY2r)]  k„{Yxr)-Y2r^niY2r)- — - 


In  solving  for  the  poles  of  the  scattered  waves,  a„  will  correspond  to  the  H-modes,  and  b„  will 
correspond  to  the  E-modes.  Therefore,  the  exact  poles,  ,  of  the  H-modes  for  a  dielectric  sphere 

can  be  found  from  the  equation 

■^[r2"«(r2'')]  K(y\^)-Y 2". (r 2'') ■  =  0 

and  the  exact  poles  of  the  E-modes  for  a  dielectric  sphere  can  be  found  fi-om  the  equation 

r  .  ,  o',  ,  .  .  ,  AYAniyx^)]  _n 

[r2"»(X2'')]  K(yxr)-Y2rin(y2r)- — ^7^: - 0 

Slab: 


We  will  now  consider  the  case  of  the  infinite  dielectric  slab,  with  thickness  / ,  in  a  dielectric 
medium.  Technically  there  would  not  be  natural  fi-equencies  since  the  pole  locations  depend  upon 
the  direction  of  the  incident  wave  [Baum,  1994].  However,  since  we  are  really  interested  in  finite 
size  targets  we  will  assume  normal  incidence.  Making  this  assumption  the  scattering  coefficient 
Rg,  has  been  derived  as  [Baum,  1994]: 

^ _ [^'-l]sinh(;'2/) 

2^cosh(;'  2O + + 1]  sinh(/  2O 
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Note  that  there  is  only  one  scattering  coefiBcient  here.  This  is  because  the  magnetic  field  is 
polarized  (parallel  to  the  electric  field,  so  the  E  and  H  poles  are  the  same).  This  implies  that  the 
exact  poles,  ,  for  the  dielectric  slab  can  be  found  from  the  equation 

2^cosh(y20  +  ['5^  +  ljsinh(f  2O  =  ^ 

Cylinder: 

For  the  dielectric  infinite  cylinder  the  we  do  not  have  the  spherical  wave-forms  that  we  had  in 
the  case  of  the  sphere.  Therefore,  the  Bessel  functions  in  the  scattering  coefficients  will  be  the 
standard  first  and  second  order  Bessel  functions  (not  spherical  or  modified  spherical).  These 
Bessel  functions  in  the  integral  form  are 

J„  (x)  =  —  J  cos(x  sin  0  -n6)d0 

^  0 

and 

H„ix)  =  J„ix)-iY„(x) 
where 

Y„{x)  =  —^s\n{xs\n9-n6)d9-—^^”*  +(-l)"e""^je"**'^'‘t/^ 

The  scattering  coefficients  a„ ,  and  derived  in  the  o  form  (s  =  jco)  are  [Balanis,  1989] 

a  ^  (y\r)J„ir2r) 

^„{Y2r)HAr,r)-JAr,r)H„{r,r) 


Jy  ^  iy2r)Jn{rxr)-J„ 

Jn  (r  2r)H„  iy  /)  -  ir  2^)H„  (r  /) 

In  the  case  of  the  dielectric  cylinder  the  a„  coefficient  corresponds  to  the  E-modes  and  the  Z>„ 
coefficient  corresponds  to  the  H-modes.  Therefore,  the  exact  poles,  £0„  (5„  ^  ©„  ),  for  the 
infinite  dielectric  cylinder  can  be  found  by 

^niy2r)Hl(yy)-JAy^r)H„{y,r)  =  0 

for  the  E-modes  and 

Jniy2r)HAY,r)-^Ay^r)H„{y,r)-0 

for  the  H-modes. 


IV.  Perturbation  Formulas 

To  the  casual  observer  these  exact  formulas  probably  don’t  look  so  bad.  However,  anyone 
who  has  ever  had  to  expand  Bessel  functions  will  tell  you  that  finding  the  poles  of  functions  that 
contain  multiples  and  derivatives  (and  multiples  of  derivatives)  of  Bessel  functions  would  be  a 
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cumbersome  task.  Recall  also  that  these  exact  formulas  are  for  nice  geometrically  shaped  objects. 
Pole  functions  for  real  ordnance  cases  would  be  even  more  busy  (much  more).  With  this  in  imnd 
Carl  E.  Baum  made  some  approximations  that  would  make  these  functions  easier  to  work  with  for 
the  slab  and  the  sphere.  The  hope  is  that  these  perturbation  functions  wiU  give  results  so  close  to 
the  exact  functions  that  others  will  be  able  to  apply  similar  approximations  to  the  exact  pole 
functions  for  the  real  ordnance  cases.  The  following  paragraphs  will  show  Baum’s  results  for  the 
slab  and  the  sphere  from  his  July  24,  1994  interaction  note  (IN  504). 


Slab: 


For  the  slab  we  will  define  the  natural  frequencies  as 

=  ^a.O  + 

where  □  ‘s  found  from  the  limiting  case  when  ^  qo  giving 

sinh(5„oV/^o^2  0  =  0 

note  that  this  number  will  be  imaginary.  For  I  refer  you  to  Baum’s  work  where 

I 

2s^2 


As  «  — 


Sphere: 


The  same  method  applied  to  the  sphere  has  similar  results  where  the  poles  are  still 
Now,  however,  we  wUl  have  different  ^ ;  one  for  the  H-modes  and  one  for  the  E-modes.  For 


the  H-modes  as  |  oo  we  see  that  the  exact  formulas  can  reduce  to 

/„  (y  jf)  =  0 ,  or  in  other  words:  /„  (s^  o  =  0 , 

For  the  E-modes  as  ^  we  can  see  that  the  exact  formulas  can  reduce  to 

\r  (r  2^)]  =  0  ’  o*’  other  words:  q  (^a.o  ^ 

Again,  o  ^  purely  imaginary  number.  The  As^  term  for  the  sphere  is  nearly  identical  to  the 
slab  where 


As 


I 


i_Ei_ 

2  ^a.O^l 


for  the  H-modes  and 


er2 


0^2^ 


1  + 


/•  0-, 

i^a.0^|Mo^2ry 

to 

o 

t  _ 

(Sa,oylMoe2f’f  +n(n  +  l) 

for  the  E-modes. 
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Cylinder: 


A  fellow  SRP  worker  (for  RDL)  George  Hanson  used  Carl’s  same  approach  for  an  infinite 
dielectric  cylinder.  The  results  again,  were  very  similar.  Where  ^  for  the  E-modes  ^ 

can  be  found  from  =  0,  or  in  other  words:  J „(,o) =  0,  and  for  the  H-modes 

0  can  be  found  from  J„  (/ 2^)  =  0 ,  or  in  other  words:  J„  (a) „  0 )  =  0 .  For  As„  George 
found 


I  r 


-e,2 


1  + 


1  O’] 

2 


for  the  E-modes,  and 


As„ 


L-Ej- 

2 


for  the  H-modes. 


V.  The  Results 


The  results  were  found  with  great  assistance  from  the  Mathcad  6.0  Plus  software  using  the 
root  finder  function.  However,  in  order  to  find  the  exact  (and  perturbed)  function’s  poles  for  the 
first  three  E-modes  and  H-modes  the  functions  had  to  be  ^eatly  simplified  before  they  could  be 
put  into  Mathcad.  This  meant  analytically  expanding,  multiplying,  and  taking  the  derivatives  of  the 
first  three  Bessel,  Hankie,  modified  spherical  Bessel,  and  modified  spherical  Hankie  functions 
(then  reducing  them).  Reducing  them  was  a  tedious  but  necessary  task  (Mathcad’s  root  finder 
wouldn’t  work  with  the  functions  in  their  expanded  form).  The  root  finder  itself  uses  the  secant 
method  (similar  to  the  Newton  method).  It  was  still  necessary  to  make  some  adjustments  to  the 
software  before  it  would  work,  and  anyone  who  has  ever  used  this  root  finder  knows  that  it  is 
highly  dependent  on  initial  guesses. 

In  the  case  of  the  cylinder  I  was  unable  to  find  the  “exact”  poles.  The  problem  revolved  around 
the  standard  Newman  functions  (and  their  derivatives).  The  software  could  calculate  values  for 
the  exact  functions,  but  to  set  them  equal  to  zero  to  find  the  roots  would  cause  the  software  to 
lock  up.  I  had  to  truncate  part  of  the  integration  involved  in  the  Newman  functions  in  order  to  get 
results.  Therefore,  what  poles  were  found  for  the  dielectric  cylinder  were  not  “exact.”  These  poles 
did,  however,  provide  a  “check”  for  the  poles  from  the  perturbation  formulas.  Note  that  I  could 
find  poles  from  the  perturbation  formulas  for  the  cylinder  but  not  from  the  exact  formulas  (hence, 
the  perturbation  formulas  are  already  proving  their  usefulness). 

In  order  to  show  all  of  my  results  for  this  project  I  would  need  the  paper  equivalent  to  that 
used  in  making  War  and  Peace.  Therefore,  I  will  instead  list  all  of  the  cases  that  I  found  poles  for 
and  present  a  graph  that  shows  the  most  effective  regions  for  the  perturbed  functions. 
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The  following  is  a  list  of  the  42  different  cases  that  poles  were  found  for  both  the  exact  and  the 
perturbed  functions  for  the  dielectric  slab  and  sphere. 
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1.5E+09 
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10cm 

1E+09 
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10cm 

Condensing  all  42  different  cases  into  a  graph  that  would  establish  an  effective  region  for  the 
perturbation  functions  required  two  jobs.  First  I  found  the  percent  errors  for  all  42  cases  between 
the  exact  versus  perturbed  poles.  Then  I  established  two  unit-less  parameters,  t  and  s^,  for 

'■O'l  ,  ^2 

graphing.  Where:  r  = -  and  e^=  — 

In  the  following  graphs  everything  inside  the  curves  represents  a  “good”  region.  In  the  first  graph 
the  poles  fi-om  the  perturbed  functions  were  less  than  6%  off  for  the  region  inside  the  curve,  and 
in  the  second  graph  (next  page)  the  poles  fi'om  the  perturbed  functions  were  less  than  2%  off 


Graph  1 :  <  6%  error  for  perturbed  functions 
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Graph  2:  <  2%  error  for  perturbed  functions 


I*  J 

Graph  3:  An  overlay  of  the  2%  and  6%  error  graphs 


VL  Discussion 

The  first  point  that  I  would  like  to  make  is  about  the  nature  of  the  poles  themselves  (exact  or 
perturbed).  The  closer  gets  to  S2  the  more  shift  to  the  left  for  the  real  components  (with  very 
little  change  in  the  imaginary  terms).  I  should  also  point  out  that  for  the  exact  poles  there  was  a 
non  complex  pole  that  doesn’t  show  up  using  the  perturbed  functions.  However,  since  real  roots 
don’t  oscillate  they  wouldn’t  show  up  on  radar  anyway.  I  did  notice,  as  theory  dictates,  a 
“layering”  effect  for  the  poles  of  the  sphere  and  slab. 
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The  effect  of  the  conductivity  term,  <7, ,  is 


j®n  ^ 

000 


O 

O 

O 

*■ 

o 

o 

o 


Sphere/Slab  Poles  for  a,  =0 


Sphere/Slab  Poles  for  increasing  <7, 


The  greater  the  conductivity  of  the  soil  the  more  acutely  angled  shaped  the  pole-layer  became. 
Where  the  perturbation  formulas  usually  had  greatest  errors  for  the  most  damped  of  those  poles 
(poles  that  were  furthest  left). 

As  for  the  size  of  the  target  usually  being  lOcm,  generating  poles  for  other  size  targets  would 
simply  involve  the  following:  if  the  target  is  increased  in  size  by  some  factor,  k,  then  the  poles 
would  be  decreased  by  this  same  factor  and  vice-versa  (for  cr,  =0;  this  is  only  an  approximation 

for  other  <7, ). 


For  cases  where  s^ »  increasing  ^ ,  by  a  factor  of  100  decreased  the  real  part  of  the  poles 
by  a  factor  of  10  (for  52=1).  The  imaginary  term  remained  the  same. 


There  were  some  unexpected  results  in  comparing  the  poles  for  the  exact  versus  perturbed 
fijnctions.  Notice  from  the  graphs  (in  the  result  section)  that  there  are  many  cases  where  some  er 
would  contain  “good”  values  but  for  (7,=0  (t=0)  the  percent  errors  were  to  large  to  be 
considered  “good.”  I  noticed  when  compiling  the  percent  errors  that  as  <7,  is  increased  (hence,  as 

T  is  increased)  the  error  decreased  until  it  reached  some  minimum,  then  the  error  would  begin  to 
increase  until  it  was  too  large  to  be  considered  “good.”  This  explains  the  reason  for  the  similar 
shapes  of  the  two  graphs  (“<6%  error,”  versus  “<2%  error”). 

The  effectiveness  region  established  for  the  perturbation  functions  would  be  applicable  to 
many  different  soils.  Where  for  soils  (with  plastic  targets)  t  ranges  between  0  and  1  and  ranges 
between  0. 1  and  0.6.  The  “<6%  error”  graph  shows  “good”  results  for  x  ranging  (nonlinearally) 
fi-om  0  to  1.75  and  ranging  (nonlinearally)  from  0  to  0.4.  The  best  results  fi-om  the  perturbed 
fiinctions  correspond  to  low-moisture  soils. 

In  comparing  the  poles  from  the  perturbed  formulas  for  the  cylinder  with  the  poles  from  the 
“near”  exact  formulas  I  was  able  to  get  low  percent  errors.  In  fact  the  results  seemed  to  show 
basically  the  same  effectiveness  region  that  was  found  for  the  sphere  and  slab  case.  However, 
since  I  was  unable  to  get  the  “exact”  poles  (due  to  the  computer  limitations  discussed  earlier)  I 
will  not  be  presenting  these  results. 
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Vn.  Conclusions 


The  exact  approach  for  finding  the  natural  signatures  for  a  dielectric  sphere,  and  an  infinite 
dielectric  slab  and  cylinder  were  presented.  This  approach  was  used  to  find  these  poles  for  various 
size  targets  and  different  dielectric  mediums. 

An  approximate  approach  for  finding  these  same  results  was  presented.  The  results  found  from 
these  approximation  formulas  were  compared  to  the  exact  poles.  A  region  of  effectiveness  was 
established  for  these  approximating  formulas.  This  region  varied  with  for  x  ranging  (nonlinearally) 
fi’om  0  to  1.75  and  ranging  (nonlinearally)  from  0  to  0.4.  The  effective  region  theoretically 

contains  many  different  types  of  soils.  The  approximation  formulas  showed  the  best  results  for 
low-moisture  soils. 

The  approximation  formulas  proved  their  usefulness  when  the  “exact”  poles  could  not  be 
found  for  a  dielectric  infinite  cylinder.  Although,  “nearly”  exact  poles  for  the  cylinder  showed  that 
the  effectiveness  region  is  close  to  that  of  the  dielectric  slab  and  sphere. 

With  these  three  very  different  shapes  having  basically  the  same  effectiveness  region,  I  feel 
fairly  confident  that  this  effectiveness  region  will  be  the  same  for  any  of  the  specific  buried  plastic 
ordnance  cases.  Therefore,  the  approximation  approach  made  by  Carl  E.  Baum  will  be  of  great 
assistance  in  developing  a  library  of  the  natural  fi'equencies  of  different  plastic  ordnance. 
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